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Abstract 



Ground state wave function of two-dimensional t — J model is found at 



I. INTRODUCTION 



Oh- 

<D ■ doping close to half filling. It is shown that the condensation of Cooper pairs 

(superconducting pairing of mobile holes) and the condensation of spin-waves 
into spin-liquid state are closely connected. The effective spectrum of S = 1 
excitations, spin- wave gap pseudo-gap Am, magnetic correlation length 
and static magnetic formfactor Sm (q) are calculated. 
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In the recent papers Jl]|| in the framework of t — J model we have demonstrated very 
strong d-wave pairing between dressed quasiholes induced by spin-wave exchange. The pair- 
ing gives critical temperature in a reasonable agreement with experimental data for copper 
oxide superconductors. In general terms our approach follows the studies, of Monthoux, 
Pines, Scarlapino, Schriffer, and others, supporting idea of the magnetic fluctuations mech- 
anism of pairing [|p§]. For calculations we have used BCS-like approximation for dressed 
quasiholes. A similar approach to the pairing of dressed quasiholes has been used by Dagotto, 
Nazarenko, and Moreo M. The important difference is that in our work the hole- hole 
interaction is derived from the parameters of the t-J model, while in || it is introduced ad 
hoc with magnitude adjusted to fit experimental data. 

We consider the holes as well as spin-waves moving in the background with long-range 
antiferromagnetic (AF) order. It means that we assume that the AF order is preserved at 
distances smaller than inverse Fermi momentum r < 1/pf- It is widely believed that this 
assumption is correct and present work confirms this point. It is also widely believed that 
correct spin structure of the ground state should be of the type of Anderson spin-liquid 
state HlQf . For t — J model the instability of long-range AF order together with normal-hole 



Fermi liquid was probably first pointed out by Shraiman and Siggia fllT]. Later it has been 



demonstrated in numerous works (e.g. see Refs. |I2]-|16|). It is due to the strong interaction 



of spin-waves with mobile holes. However, a structure of the ground state as well as a 
spectrum of excitations was not understood, and spin-liquid state was not derived. 

Now we understand that we have to consider a hole-Fermi-liquid with superconducting 
pairing which interacts with long-range fluctuations of AF ordering. Physical difference of 
the paired Fermi-liquid from normal one is stiffness of the former. In the present work we 
demonstrate that the long-range AF order is also destroyed, but pairing is preserved. We 



derive explicitly the spin-liquid state of the system. We show that selfconsistent translation 
invariant spin-liquid solution exists at doping 5 > 5 C ~ 3 — 4%. Using this solution we 
calculate spin- wave pseudo-gap Am, magnetic correlation length £m an d static magnetic 
formfactor 5(q). In general terms the results which we derive here in the framework oft — J 
model are very close to the picture of magnetic state discussed in the works of Chubukov, 
Pines, Sachdev, Sokol, and Ye P?|-pO[] . 

The paper has the following structure. In Sec. II we present the effective Hamiltonian of 
t — J model and discuss pairing of holes. In Sec. Ill we formulate the idea of our approach pre- 
senting the simplified calculation of spin- wave Green's function. Sec. IV presents calculation 
of spin-wave polarization operator. In Sec.V spin-wave Green's function is calculated. We 
also discuss the pseudo-gap in spin-wave spectrum. Sec. VI presents calculations of magnetic 
correlation length £ M and static magnetic formfactor ^/(q). Finally our conclusions are 
given in Sec. VII. 



II. EFFECTIVE HAMILTONIAN AND PAIRING OF DRESSED HOLES 

The t-J model is defined by the Hamiltonian 

H = H t + Hj = -t + H.c.) + J J2 Sn-Sm, (1) 



<nm>(T <nm> 



where d\ a is the creation operator of a hole with spin o io =t? I) at site n on a two- 
dimensional square lattice. The d\ a operators act in the Hilbert space with no double electron 
occupancy. The spin operator is S n = ^S na cr a pd n p. < nm > are the nearest-neighbour sites 
on the lattice. Below we set J as well as lattice spacing equal to unity. 

At half- filling (one hole per site) t-J model is equivalent to the Heisenberg AF model 
2I]fZ3] which has long-range AF order in the ground state [E3,E3]. Let us denote the wave 



function of this ground state by |0). This is undoped system. We consider the doped system 
basing on the ground state of undoped one. In spite of destruction of the long-range AF 
order it is convenient to use |0) and corresponding quasiparticle excitations as a basis set 
in the problem with doping. The effective Hamiltonian for t — J model in terms of these 
quasiparticles was derived in the papers 



H eff = Yl e k^L^k CT + ^q( a q«q + 0qAl) + H h,sw + H hh . (2) 



kcr 



It is expressed in terms of usual spin-waves on AF background a q , /3 q (see e.g. review P% 



and composite hole operators (er = ±1/2). The summations over k and q are restricted 
inside the Brillouin zone of one sublattice where 7 q = |(cosg x + cosg y ) > 0. Spin- wave 
dispersion is 



cu q = 2JI - 72 Wq «V2|q|, atq<^\. (3) 



Single hole properties are well established (for a review see Ref. [29fl). Wave function of 



a single hole can be represented as ipk a = h^O). At large t the composite hole operator 
has complex structure. For example at t/J = 3 the weight of bare hole in ip^ is about 
25%, the weight of configurations "bare hole + 1 magnon" is ~50%, and of configurations 
"bare hole + 2 or more magnons" ~25%. Dressed hole is a normal fermion. Hole energy 
ek has minima at k = k , where k = (±7r/2, ±7r/2). For t < 5 the dispersion can be well 



approximated by the expression [30 



e k ss E + 2 - ^0.66 2 + 4.56t 2 - 2.8t 2 ~fl + ^ 2 (cos k x - cos k y ) 2 . (4) 

The decimals in this formula are some combinations of the Heisenberg model spin correlators. 
Constant E defines reference level for the energy. To find e k with respect to undoped system 
one has to set Eq = 0. However for present work it is convenient to set ek = 0, and therefore 
Eq = V0.66 2 + 4.56i 2 — 2. The coefficient ft is small and therefore the dispersion is almost 
degenerate along the face of the magnetic Brillouin zone 7k = 0. According to Refs. [31 ,[32 



ft ~ 0.1 • t at t > 0.33. To avoid misunderstanding we note that formula (|J) certainly is 
not valid for very large t where the hole band width is saturated at the value of the order 
of unity and does not increase with t. However physically we are interested in t « 3 (e.g., 
see Refs. j33H35j) where @) works well. Near the band minima k the dispersion @) can be 



presented in the usual quadratic form 

e P ~ \PiPl + 2&#L P2 < ft, (5) 

where pi (P2) is the projection of k — k on the direction orthogonal (parallel) to the face of 
the magnetic Brillouin zone (Fig.l). From Eq.(|) for t > 0.33 we find ft « 0.65t, hence the 
mass anisotropy is ft /ft ~ 7. 

The interaction of a composite hole with spin waves is of the form (see, e.g., Refs. 



H h>sw = #k,q ( ^k+ql^kTttq + 4+qT^I Al + H.C. J , (6) 

k,q V ' 

# k ,q = 2v / 2/(7kf/q + Tk+q^q), 



where Uq = J-^ + \ and Vq = —sign(jq)^^ — \ are the parameters of Bogoliubov trans- 
formation diagonalizing spin-wave Hamiltonian. The hole spin-wave coupling constant / is 
a function of t evaluated in the work [p5 ]. For large t the coupling constant is t-independent 
/ ~ 2. Let us stress that even for t > J the quasihole-spin-wave interaction (||) has the same 
form as for t ^ J (i.e. as for bare hole operators) with an added renormalization factor (of 
the order J/t for t ^> J). This remarkable property of the t-J model is due to the absence 
of single loop correction to the hole-spin-wave vertex. It was first stated implicitly by Kane, 
Lee, and Read P7|. In Refs. pl|,p6|,p5| it was explicitly demonstrated that the vertex cor- 



rections with different kinematic structure are of the order of few percent at tj J ~ 3. There 
is also some g-dependence of the coupling constant /. For example f(q = it) w 1.15 f(q = 0) 
at t/J = 3 (see Refs. p7|j2]|). However this dependence is weak, it is practically beyond 
the accuracy of the calculation of the renormalized value of /. Therefore we neglect this 
dependence. 

Finally there is contact hole-hole interaction in the effective Hamiltonian (|2|). is 
discussed in details in the Refs. [p6 , |27 ,2|. It is proportional to some function A(t). For small 



t this function approaches the value -0.25, which gives the well known hole-hole attraction 
induced by the reduction of the number of missing antiferromagnetic links. However for 



realistic superconductors t ~ 3 (see e.g. Refs. |33H35|). Surprisingly function A(t) vanishes 



exactly at t ~ 3, and it means that the mechanism of contact hole-hole attraction is switched 
off. In contrast the spin-wave exchange mechanism Hh, S w is negligible for small t where / ~ t, 
and it is most important at large t where / approaches 2. We are interested in "physical" 
values oft: t ~ 3. Therefore in the present work we neglect contact interaction (H^ = 0) and 
consider only the hole-spin-wave interaction at t — 3. Corresponding value of / according 



to @] is / = 1.8. 



The interaction which gives hole-hole superconducting pairing arises from spin-wave ex- 
change (Fig.2). From (Q) we find 

if o fl'k,qfl , k',-q /~\ 

V W = ~ l — p p~> \'l 

—LUq — £/ k — _C/ k / 

where q = k + k', and i? k is the energy of quasihole excitation above the BCS ground 
state (see below). We consider the interaction in static approximation. Justification of this 
approximation is as follows. At small hole concentration 5^1, the holes are localized in 
momentum space in the vicinity of the minima of the band ko = (±7r/2, ±7r/2) and Fermi 
surface consists of ellipses. The Fermi energy and Fermi momentum of non-interacting holes 
are 

e F « ^{PifafH, p F « y/p^ « (vr^) 1 / 2 . (8) 

The Fermi momentum p F is measured from the center of corresponding ellipse. In pairing, 
the exchange of spin- waves with typical momenta it > q > p F is the most important (see 
also Refs. |]J|,|27j ) . The energy of such spin wave is much higher than the typical energy of 
a Cooper pair 

uJg ~ q > p F ~ (vr5) 1/2 > e F ~ (/3i/3 2 ) 1/2 5. (9) 

This justifies the static approximation (|7|) for hole-hole interaction (see also comment at 
end of section V). This is much different from usual phonon induced pairing where Debye's 
frequency is much lower than the Fermi energy. 
With trial function of ground state 

w=nK+^4 T ^ki)io>. ( io ) 

k 

we get conventional BCS equation for superconducting gap A k 

A k = -iE^^tanh^, (11) 



£k 

M k 




Vk = signAi 

where ^ k = e k — /x, and \i is chemical potential fixed by equation 



*=2E(«4 + £^t) (12) 



All details of approximate analytical and exact numerical solutions of equation (|TT| ) are 
discussed in our papers This equation has an infinite set of solutions of different sym- 
metries with strongest pairing in the d-wave (four nodes of the gap at the face of magnetic 
Brillouin zone). For further discussion we need only the results for d-wave. All numerical 
values are calculated at t = 3, P1/P2 = 7, and / = 1.8. In Table I we present for different 
values of hole concentration l)Fermi energy e F , 2)chemical potential at zero temperature 
/i(0), 3)maximum value of the gap in Brillouin zone at zero temperature A max (0), ^maxi- 
mum value of the gap at Fermi surface at zero temperature A Fmax (0), 5)chemical potential 
at critical temperature /i(T c ), 6)Critical temperature T c . In Fig.l we give the map of a single 
hole mean occupation number f 2 . at 5 = 0.1 and zero temperature. 



III. MODIFIED SPIN-WAVE THEORY AND QUANTUM MELTING OF 
LONG-RANGE ANTIFERROMAGNETIC ORDER 



It is well known that in the long wave-length limit Heisenberg model is equivalent to 
the nonlinear cr-model (e.g., see review paper f28fl). Therefore usual field theory crossing 
symmetry is valid. The same is valid for t — J model. Technically it is evident from Eq.(|6]): 
at q <C 1 the vertex g k q pa <?k-q,q- This means that instead of considering a set of Green's 
functions (a^al), (ctk/3-k)---, one can introduce one combined Green's function of vector 
excitation 

^ ' ^ oj 2 — u 2 — 2u; q n(u;, q) ' ^ 
where H(uj, q) is mobile holes polarization operator. For stability of the system the condition 

cu g + 211(0, q)>0 (14) 
should be fulfilled. Otherwise the Green's function (O) would possess poles with imaginary 



to. It was demonstrated in the papers piHlql) that with 11(0, q) calculated in the normal hole 
Fermi liquid approximation the condition (|TJ|) is violated. It means the instability of long- 
range AF order. Now we want to: l)Take into account strong hole-hole pairing, 2)Formulate 
the approach for description of state without long-range AF order. In the present section 
we discuss only point 2). 

For elucidation of physical meaning of our approach it is convenient to use Hamiltonian 
technique instead of the Feynman one. Due to the interaction with mobile holes the wave 
function of the renormalized spin-wave corresponding to the Green's function fll3|) is a com- 
bination of a q and /3_ q . To find this wave function let us write down the effective spin-wave 
Hamiltonian. 

H s W = W(u; q + n(u;,q))(a q a q + /3 q /3 q ) - n(u;, q)(a q /3_ q + a q /3l q )Y (15) 

q \ ' 

The term proportional to co> q comes from the "bare" Hamiltonian (fj). First "n-term" comes 
from the diagram Fig. 3a where one spin-wave is annihilated and the other is created. Second 
"n-term" comes from diagrams Fig.3bc where two spin-waves are annihilated or created. Let 
us note that spin- waves have definite values of S z : a q has S z = — 1 and /3_ q has S z = +1. 
Therefore they can appear only in combinations presented in flT5|). One can certainly prove 
this explicitly using (||) and calculating the polarization operator. In the second "n-term" 
the spin-waves have the opposite momenta. At q <C 1 the vertex is proportional to the 
momentum q. Just for this reason second "n-term" has different sign in comparison with 
first one. Diagonalization of the Hamiltonian (|i~5f ) by the Bogoliubov transformation gives 
the spectrum of Bose excitations in the system 

fi q = u4 + 2u; q n(fi q ,q). (16) 



This is exactly the equation for the poles of the Green's function fll3|). So the Hamiltonian ap- 
proach reproduces conventional result of the Feynman technique. Account of tu-dependence 
of the polarization operator in Hamiltonian technique is rather cumbersome. Therefore in 
the present section we neglect it: n(fi q , q) pa n(0,q). As usual Bogoliubov transformation 
gives the ground state of the form 

|^)ocexpfec q a q /3O|0), (17) 



with some coefficients c q . This is the condensate of spin- waves. 

We started from Neel ground state |0) with two sublattices w-up and <i-down. The 
difference in magnetization of the two sublattices is of the form (see e.g. Ref. p§| ) 



\{S Z U - S z d ) = 0.303 - 2 Y, — («i«q + $A - 7 q (« q /3- q + atyl 



q 



q 



Using parameters of transformation which diagonalizes the Hamiltonian (|15|) one can easily 
calculate the renormalized magnetization 

18. = - Sa M = 0.303 - 2/(J- - J-)|S_, (19) 



where f2 q = y ^ q + 2cj q n(0, q). 

In the above discussion we assumed that condition fll4|) is fulfilled and hence f2 q is real. 
Let us now increase polarization operator (11(0, q) — > x ■ n(0, q), x > 1) approaching O q to 
zero. There are two possibilities: 1) f2 q vanishes at some values of q, but 5S Z remains finite 
positive because of the convergence of the integral in (|T^). 2) In approaching f2 q to zero 



5S Z vanishes and then becomes negative. The choice between these two scenario depends 
on q-dependence of 11(0, q). We will demonstrate below that for paired hole Fermi liquid at 
5 > 5 C second scenario is realized. Therefore let us discuss this situation. 

Vanishing of the magnetization 5S Z means that there is a lot of spin-waves in the con- 
densate ([[7]) and we have to take into account their nonlinear interaction. We can not do 
it exactly. Fortunately, there is approximate way. We can apply the modified spin-wave 



theory suggested by Takahashi for the Heisenberg model at nonzero temperature fl38fl . Fol- 
lowing Takahashi let us impose the condition that the sublattice magnetization vanishes in 
the quantum spin-liquid state 

(9s\l(S z n -S z d )\gs) = 0. (20) 
To find the ground state with this condition we have to diagonalize 

H v = H sw -\v 2 {S z u -Sl), (21) 



where %u is Lagrange multiplier, H sw is given by (|l5l) and (S* — S%) by (|T|). Simple 
calculation shows that instead of (|16|) we get a spectrum of excitations with a gap 



^q = V fi q + y2 - ( 22 ) 

The average value of the magnetization is given by the formula (|19|) with f2 q replaced by 
f^q. The gap v should be found after substitution of this formula into condition (|20|). Let us 
stress that this condition reflects strong nonlinearity of the spin-wave theory. In essence it 
gives an effective cutoff of unphysical states in the Dyson-Maleev approach. This question is 
discussed in the paper |39| where the modified spin-wave theory is applied to the description 
of the transition from AF state to spin-liquid state in J\ — J 2 model. 

In the conclusion of present section we would like to note that realization of the first 
scenario (finite positive magnetization 5S Z at vanishing of O q at some values of q) would 
mean instability with respect to decay to spiral state. It happens at 5 < 5 C . 



IV. SPIN- WAVE GREEN'S FUNCTIONS AND POLARIZATION OPERATORS 



We are interested in hole concentration 5 ~ 0.15. Corresponding value of q is q ~ pp ~ 1. 
Therefore, for practical calculations cx-model long-wave-length approximation ([13]) is not 
enough. The technique which overcomes this problem was developed by Igarashi and Fulde 
40|| , see also the paper by Khaliullin and Horsch El] . Following these works let us introduce 



the Green's functions for spin- waves 

D ao (t > q) = -i(T[a q (t)at(0)]), (23) 
D a p(t,q) = -i(T[a^t)/3^{0)]), 
D/9a(*,q) = -i(T[/3i q (t)at(0)]), 
D^(t,q) = -i(T[/3i q (t)/3_ q (0)]), 

where T is the time-ordering operator, and (...) represents an average over the exact ground 
state. In zero approximation in the interaction with mobile holes the spin-wave Green's 
functions in u-q representation are given by El 



D° aa (u,q) = (uj-u J(l + zO)-\ (24) 
D° ap (u,q) = D° pa (u,q) = 0, 
D%(uj,q) = (-lj-^ + iO)- 1 . 

Polarization operator in this technique is also a matrix. So there are P aa (uj,q), P a p(u>,q), 
P/3 a (u>,q), and Ppp{uj,q) components. 

In one loop approximation polarization operator is given by the diagrams presented in 
Fig. 4. where the normal as well as the anomalous hole Green's function is taken into account. 
First of all let us discuss contributions of the "coherent" and "incoherent" parts of the hole 
Green's function into the polarization operator. The diagrams in Fig. 4 are represented in 
terms of the dressed quasiholes which involve in the effective Hamiltonian (Q). However 
one can ask the question: what is the accuracy of this approximation and what happens if 
we substitute the Green's function of bare hole d na into polarization operator? The Green's 
function of bare hole is of the form 

G(e, k) = _ Z + G mcoh . (25) 

where Z is quasiparticle residue. If we use bare hole Green's function we also have to use 
bare hole-spin- wave coupling constant f bare = It in the vertex g kq (^|). Substitution of the 
pole (coherent) part of (|25|) into polarization operator gives combination fbareZ, but this is 
exactly the effective coupling constant / ||25|| . Thus, the effective theory with Hamiltonian (0) 



is equivalent to the account of the coherent (pole) part of the bare hole Green's function ( 25f) . 
According to condition (113) the most important characteristic is the polarization operator 
at zero frequency. If we neglect the pairing of quasiholes we can easily calculate its value in 
normal Fermi-liquid approximation fi~6|| ) 

/2 f 2 

P aQ (0,q)=n(0,q)^-^=g, at q<^p F . (26) 

7TVPlP2 

This is coherent contribution. It is independent of the hole concentration! Incoherent part 
of the polarization operator which comes from Gi nco h in (f25|)is proportional to hole concen- 
tration 5. Therefore it is negligible at 5 <C 1. This conclusion agrees with that of the Ref. 



42j. It is interesting that incoherent contribution can be also calculated analytically 



P t ncoh(0,q)^--^5\n]-q. (27) 

At any reasonable 5 the incoherent part Pi nC oh(0,q) does not exceed 10% of P co / l (0,q). 
Therefore further we neglect Pincoh and use effective theory with the Hamiltonian (g) which 
is equivalent to the account of only coherent part in the hole Green's function (|25|). 

Due to the symmetry of Green's function it is convenient to introduce following notations 

Ppp{-u, q) = P aa (uj, q) = U(u, q), (28) 
P a p(u, q) = Pp a {v, q) = n(w, q). 

Calculation of the polarization operator Fig.4 in single loop approximation is straightforward. 
Account of both normal and anomalous hole Green's functions gives 



n (1 V,q) = E 



,2 „2 



k L \UJ — -C/ k — -C/k+q — ^ — -C'k — ^k+q/ 

/ 1 1 

-ttkVktt-k-qV-k-q^k,qPk+q,-q p p 1" 



X 



^ — E k — -^k+q — U) — E k — E k+ q 

1 '' ' U'. qj = 

k 

f 1 1 

x ^ ^ + 



n (1 V 5 q) = XA U k W k+q#k,q#k+q,-q ~ UkVk«-k-qU-k-q0k, q 



\0; — — -Ek+q — ^ — -^k — -^k+q 

In these formulas (?k, q is hole-spin- wave vertex @, and «k, fk, -Ek are parameters of BCS 
hole wave function (p"Tp. 

Let us consider now two and more loop corrections to polarization operator. Some 
diagrams of this type with normal hole Green's function are presented in Fig. 5. One can 
easily prove that the diagram Fig. 5a (single spin-wave exchange) vanishes due to the spin-flip 
nature of the hole-spin-wave vertex. Vanishing of this diagram is absolutely similar to the 
vanishing of the single loop correction to the hole-spin- wave vertex (see Refs. [31, 36. 25]). 



The diagram Fig.5b equals to zero for the same reason. To calculate the diagram Fig. 5c let 
us represent it using Goldstone diagram technique: Fig. 6. In the intermediate state in the 
diagram Fig. 6a as well as in the diagram Fig. 6b there is two-particles two-holes excitation 
HJ. Nevertheless these two diagrams are essentially different. All intermediate momenta 



in the diagram Fig. 6b are in the vicinity of Fermi surface (p = k — ko ~ Pf), and therefore 
this diagram is convergent. On other hand the particles momenta in diagram Fig. 6a are 
not restricted and this diagram is ultraviolet divergent. The divergence gives logarithmic 
enhancement of this diagram. Further we neglect contribution Fig.6b, and transform Fig. 6a 
into an effective diagram Fig.7 with "dot" given by Fig.8. 

We are interested in relatively small q (q < 1). Therefore all momenta in the diagram 
Fig.7 are close to the Fermi surface. It means that we need to calculate the effective in- 
teraction ("dot") at Fig.8 at k w k + q w (±tt/2, ±7r/2) and k' w k' + q w (±vr/2, ±tt/2). 
Let us neglect q-dependence (set q = 0) and denote this effective interaction by V^(k, k'). 

Single particle state near k = (— n/2, — n/2) is related to that near k 0++ = (7r/2,7r/2) 

via umklapp process and therefore these points are actually equivalent. It is like one pocket 
of the Fermi surface. Due to this reason o++ , k — ) . Certainly 

one can prove this relation explicitly using perturbation theory expression corresponding to 
diagram Fig.8. So there are only two independent pockets of the Fermi surface (for example 
centred near k 0++ and k 0+ _), and the effective interaction can be characterized by the two 
values only: 



VCS = V (1) (ko f+ ,ko++), (30) 

vl 1 l = v^(k 0++ ,k 0+ -)- 

First one is the effective particle-hole interaction when particle and hole are from the same 
pocket. Second value is particle-hole interaction when particle and hole are from different 
pockets. 

Using hole-spin-wave interaction (f|) one can easily derive expression for V++ correspond- 
ing to diagram Fig. 8 

2 2 (a 2 \ 

T/W _ QOf4V^ 7Q+k7Q-kl 4 - ^q) , , 

Q ^Ql^Q + -^Q+kJ -&Q+k 

Here k = k 0++ . At ]9p < Q < 7r the integral in fl31~|) is logarithmic divergent (/ dQ/Q) and 
therefore V++ can be analytically calculated with logarithmic accuracy 

^■gfl+Vgftj.t, (32) 



"ft (l + v/ft/A) 2 



The big logarithm L equals to 



L = ln(e A /A), (33) 



where A is typical value of superconducting gap and e\ ~ 2 is ultraviolet cutoff which is 
of the order of hole band width. Actually only this logarithm justifies introduction of the 
effective point-like interaction V^ l \ Certainly expression ([H]) can be easily exactly evaluated 
numerically. Similar consideration shows that due to the structure of hole- spin- wave vertex 
(^j effective particle hole interaction for particle and hole from different pockets vanishes 

Vl 1 } = (34) 

Above discussion is not the end of the story about effective particle-hole interaction in 
the spin-flip channel. It was just a leading order. To find correct value we have to sum all 
the "ladder" given at Fig. 9. Fortunately each additional rung in ladder gives additional big 
logarithm and therefore summation can be easily done with logarithmic accuracy 

fl(l + f - 



V ++ = Vmi + J • L j . (35) 

Summation of the ladder does not change the relation (0): effective particle hole interaction 
for particle and hole from different pockets vanishes. 

For very small hole concentration 5 the superconducting gap A is approaching zero, 
L — > oo, and hence 

^ _ 2/2 i±Mia (36) 

In practice L is never big enough to reach this limit. Realistic value of L can be found by 



comparing (|32|) with the result of exact numerical calculation of (|31|). At 5 = 0.05 L « 4, 
and it drops down to L w 3 at at 5 = 0.2. Substituting correct value of L into Q55D we 
find renormalized particle-hole interaction in spin-flip channel V ++ . Numerical values of V ++ 



are presented in the Table. II. Actually it is almost constant V ++ ~ —4.5 for interesting 
concentrations. 

Now we can calculate exact spin-wave polarization operator chaining effective particle 
hole interaction in the spin-flip channel, see Fig.fO. This chaining is substantially simplified 
by the fact that the pockets of the Fermi surface are well separated. Let us represent single 
loop polarization operator ( p9[ ) as a sum of contributions corresponding to different pockets 

U^ = U { P+U^\ (37) 
ke/ 

ni 1 / ) (^,q)= ^%k,q), 

kell 

where X(u,\s., q) is the integrand in formula (|29|). Regions I and II of magnetic Brillouin 

zone are shown at Fig.l. We remind that point ko = (— tt/2, — vr/2) is connected with 

k 0++ = (tt/2, tt/2) via umklapp process and therefore the region / represents one pocket. 
The same is valid for II. Due to the condition (0) particle-hole scattering does not transfer 
particle-hole excitation from one pocket to another, but inside pocket scattering amplitude 
(p5|) is constant. In this situation summation of chain in Fig. 10 is trivial. The exact spin- 
wave polarization operator fl28|) is of the form 



1 ' ^ 1 + V ++ QM q) + 1 + V ++ Q n {u, q) ' ^ 



1 + V ++ Qj(uj iq ) 1 + V ++ Q n (u,qy 
here 

Q ItII (uJ, q) = f W kWk+q + Mk^kM-k-q^-k-q) ( 7-1 1 7-1 + £T p ) • 



(39) 

Unfortunately further analytical calculation of H(uj, q) and TL(u, q) given by (^) is im- 
possible. However numerical calculation using (|29|), fl3"7|),(|3DD and parameters of supercon- 
ducting wave function found in the Section II is straightforward. We will discuss it in detail 
in the next section. Let us look here only at the limit q — > 0. In this limit 

1%, q) = U{-u, q) = -n(o;, q). (40) 

This reflects the fact that we can forget about all complications connected with multi compo- 
nent spin- wave Green's function (pH), and use simple a-model picture (|13|). Let us introduce 
parameter of spin-wave instability 

c= -2n(o,g) 

^q 



According to equation (111) if C > 1 spectrum of excitations is unstable. Naive value of C 



corresponding to normal Fermi liquid approximation (|26l) is 



C„ = -4== « 2.8. (42) 

7T\/PlP2 



Numerical values of C obtained with (|38j) for different hole concentrations are presented 
in the Table II. We see that pairing and particle-hole rescattering (|35|) substantially reduce 
value of C. Actually rescattering is more important, but it has the same origin as the pairing. 
Nevertheless C > 1, and long-range antiferro magnetic order is unstable. 



V. SPIN- WAVE GREEN'S FUNCTION. 



To deal with spin-wave instability let us apply approach developed in the Section III. 
With Takahashi condition (f20|) , (|I~8"D bare spin- wave Hamiltonian is transformed to 



q q 



E 









- " 2 ) 



z/ 2 

;« q a q + /3 q /3 q ) - 7 q - (o!q/3_ q + a q /£ q ) 

zu; q 



(43) 



q 

Lagrange multiplier v will be found later from the condition (^). Dioganalization of the 
Hamiltonian (|43"D by usual Bogoliubov transformation gives new bare spin-wave spectrum 



^4 

< = W q + ^ + -- W q + ^» ( 44 ) 



16 

and new bare spin-wave operators A' , i? q 



ttq = C q A q + Sqfilq, (45) 
/3-q = Sq4_ + C q fi_ 



q °q^q 1 l q^-q) 
Cq, >S C 



LU u(l ± Wq 



(46) 

We will see below that v is very small. Therefore we neglect i/ 4 /16 in comparison with v 2 . 
Hole-spin-wave interaction (|D should be also expressed in terms of A q and £? q operators 

Hh,sw = X/ ^ k -q ( 

q^kT^q + 4+qT^q + H.C. , (47) 
k,q V J 

Gk.q = C q (?k,q + Sq^k+q-q- 

Let us note that due to the Goldstone theorem gk, q — ► at q — > 0. However, in the state 
which we discuss now rotational symmetry is restored and therefore Gk,o = 1 /2/7k- 

Similar to (^3) we introduce Green's functions Daa^-, q), D^ B (t, q), D B A(t, q) and 
DBB(t,q). Zero approximation Green's functions D^(u;,q)... are given by eqs. (f2~4l) with 
replacement u; q — ► cov q . Polarization operators are given by the formulas of the previous sec- 
tion with replacement gk,q ~~ > ^k,q- Spin-wave Green's functions obey usual Dyson equations 



40.41 



Daa = P>aa + DaaPaaDaa + DaaPabDba, (48) 
Dbb = Dbb + D BB P BB D BB + DbbPbaDab, 
Dab = Da A PaaDab + DaaPabDbb, 
D B a = DbbPbaDaa + DbbPbbDba- 
Solution of these equations together with (|8|) gives 

D AA {^ q) = D BB (-U, q) = ~«> ~ K(-U>) (49) 

D A B{u,q.) = D BA (u,qj 



A(w,q) 
A(w,q) = 



uZ + u^+u^q n(w) + n(-w) +w n(w)-n(-w) + 1%) • n(-w) - n (w). (so) 



r2 , 



In these formulas we omit for simplicity argument q in the polarization operator. We do not 
specify imaginary part of A in the poles because below we perform Wick's rotation uo — > z£ 
which automatically gives correct behaviour. Using fll8|) one can easily express average 
magnetization in the terms of Green's functions 



= (\{S Z U - SD) = l-2i J2—\D A A(t = -0, q) + D BB (t = -0, q) - 

* q ^q 



- 7q D AB (t = -0, q) + D BA (t = -0, q) = (51) 



1 - 2 ^?f = - 303 - 2 ^(?r-i7 

q 1L L/q q ^^'l/q ^q 



where fl U(l is given by the expression 

1 1 r+°° d£ 1 / — \ 

<T = — ^UJ^S 2^q+[n(^q)+nHe,q) + 2 7q n(ie,q)] ■ (52) 

We have to find the Lagrange multiplier v from equation (|5l|) . At small q, due to the relation 
(|40|), the expression in square brackets in (Q) vanishes and A(i£, q) — > £ 2 +^ q +2coVqII(i£, q). 
If one neglects also ^-dependence of polarization operator, the integration over £ in (|5^) is 
trivial and it is reduced to formula (^) derived in Section III in a-model approximation. 

Now we can proceed to the discussion of the results of numerical calculations. Dependence 
of hole-hole pairing on the modification of spin- wave spectrum is rather weak at S <C 1, and 
we neglect it (see also end of these section). Therefore, we use parameters of hole BCS wave 
function found in Section II. Calculation of the polarization operator is described in detail 
in Section IV. One should not forget about the replacement (?k, q - ► but numerically it 
is not very important. The values of v found from equation ([|l]) are presented in the Table 
II. The effective spin-wave frequency f2j, q depends both on direction and magnitude of q. It 
is convenient to take the bare spin- wave frequency u; q as an argument instead of |q|. The 
plots of f2„ q as a function of u; q for the directions q oc (1, 0) and q oc (1, 1) and for different 
hole concentrations are given in Figs.lla-d. 

At very small hole concentration 5, equation Q5"T|) has no solution. One can easily un- 



derstand the reason for this. The parameter of spin-wave instability C is larger then 1 at 
arbitrary small hole concentration 5, see eqs.(f4lD,(f42D and Table II. It means that depen- 
dence of Q m on |q| should be of the type of solid line at Fig.lla with minimum at go ~ Pf- 
Let us denote the value of Q ucl at this point by flo. This minimum value is directly related 
to v. Therefore instead of v we can consider flo as a variable which should be found from 



eg. (|5TD - If frequency f2,,q is isotropic fi^ q ~ Q% + a(q — q ) 2 the equation ([|l]) has obvious 



solution at arbitrary small 6 

Q Q ~ e - 1/qo ~ e- 1/pF ~ e- 1 '^ (53) 

The problem is that the polarization operator and hence f2^ q is not isotropic. It is seen 
explicitly from Fig. 11a. The anisotropy is proportional to some power of go \f5. Therefore 
at small enough concentration the anisotropy is larger than exponentially small Qq and 
equation (|5l| ) has no solution. Thus if we start from spin-liquid phase (say from 5 = 0.1) 
and decrease hole density 5, at some critical value S = 5 C the gap f2 vanishes. It means 
instability of spin liquid state with respect to decay to spin waves with momentum qo 
(spirals). The critical value 5 C ~ 3 — 4%. Due to the strong exponential dependence of (|53|) 
on S it is hard to find numerically the exact value of 5 C . We would like to note that the dip 
in Fig.lla is a result of virtual admixture of spirals to the spin-liquid ground state. 



One can consider v as a gap Am in the spin- wave spectrum at q = 0. At small hole 
concentration (5 = 0.05,0.1) minimum of the effective spectrum is shifted from the point 
q = 0, see Figs.llab. One has to be careful with interpretation of Am- The effective 
frequency f2^ q defined by eq.(^) is not a frequency of simple excitation in the sense that 
Green's function G oc 1/{uj 2 — fi^ q ). Green's function is of essentially more complicated 
form. It has poles and cuts. Therefore, using standard terminology, it is better to say that 
Am is pseudo-gap and that Q m is pseudo-spectrum. The calculated value of pseudo-gap as 
well as its dependence on the hole concentration reasonably agrees with experimental data 



In equation ([51]) we have overcome the problem of real spectrum by integrating over 
imaginary to. For real u structure of the spectrum is as follow. At very small q <C pp 
there is a damped spin-wave with frequency u ~ Am ~ v- The decay of the spin-wave into 
particle-hole pair is allowed due to the nodes of the superconducting gap. With increasing of 
q, spectrum is split into the two branches: l)damped spin-wave with frequency to ~ w q , and 
2) low frequency collective excitation consisting mainly of the particles and holes. Starting 
from some value of q decay of this excitation into particle and hole is forbidden due to the 
superconducting gap. Therefore in this low frequency region spin-wave Green's function can 
be represented as 

G(u,q} = ^-^, (54) 

where Z is the spin-wave residue. The frequency of collective excitation Oq clS db function 
of co> q for direction q oc (1,0) and for hole concentration 5 = 0.1 is plotted at Fig. lib. 
Due to the decay to real particle-hole pair this spectrum has no meaning at very small q 
and at large q. Spin- wave residue Z changes along the drawn curve from Z = 0.75 at left 
hand side of curve to Z = 0.15 at right hand side of the interval. Detail analysis of the 
spectrum of spin 5=1 excitations is an important problem. It requires calculation of 
the dynamic magnetic formfactor SM(u,q) of the system. Having this formfactor one can 
perform accurate comparison with an experimental data on neutron scattering and NMR. It 
is a subject of a separate work. Here we concentrate on the magnetic correlation length £m 
and static magnetic formfactor iSjif(q). 

This is the place to remind the reader that in the calculation of the superconducting 
pairing we have used unperturbed spin-wave spectrum. This approximation is justified by the 
parameter \/5 <C 1 (see Refs. and eq. ([])). However for most interesting concentrations 
(5 ~ 0.15 — 0.2) this parameter is not so good. We see from Fig. 11a that for S = 0.05 the 
effective spin- wave spectrum is actually close to the unperturbed one. However for 5 = 0.2 
(Fig. lid) deviation is rather big. Therefore to complete our program we have to solve 
Eliashberg equations for pairing with renormalized spin-wave Green's function (compare with 
Ref. |l9||). This should be a subject of a separate work. However qualitatively the influence of 
this effect is evident. Effective pairing interaction (|7|) is roughly proportional to the inverse 
spin-wave frequency. Therefore, the renormalized interaction gives stronger pairing than 
that in "zero approximation" . Due to our preliminary estimations, the correction to pairing 
is about 30% at 5 = 0.2. 



VI. MAGNETIC CORRELATION LENGTH AND STATIC MAGNETIC 

FORMFACTOR 



We remind the reader that the spin wave operators a q and /3 q are introduced using 
Dyson-Maleev transformation for localized spin S, 



S t = a], = (2S — ajai)ai, 

Si = S — a\ai, for I e up sublattice] (55) 

S m = b m , S rn = b m (2S — bljjjm), 

Sm = ~S + HnPmi f or m £ dawn sublattice. 

The operators a q and /3 q are related to a\ and b m by usual combination of Fourier and 
Bogoliubov transformations, see e.g. Refs. |38]|28|. Following Takahashi |38| we introduce 
notations 

f(ri -i>) = (a\ai>), l^l', 
f(r m - r m /) = (blfimr), m ^ m' , (56) 
g(r t - r m ) = {a\b ] m ) = (a z 6 m ), 

where brackets (...) represents an average over ground state. Using mean field procedure 
for the averaging of the quartic terms ((a^atfb) = (a)a)($b) + (a Jf b Jf )(ab)) one can express 
spin-spin correlators in terms of functions / and g (see Ref. pg[|). 

(SjSj) = f 2 (ri — Tj) — ~<5i,j, i,j £ same sublattice, (57) 
(SjSj) = —g 2 (ri — Tj), i,j £ different sublattices. 



Using definition (|56"D, the functions / and g can be expressed in terms of Green's functions 
Daa{1 = —0, q), DAsit = — 0, q) ... This is quite similar to the average magnetization (|5T|) . 
Then with the help of formulas (|49|) one gets 



„«qr 

/( r ) = 2E7T-> 

q ii t / q 

jqr 

*00 = -2 £ Sr> 

q "vq 



The effective frequency f2 yq is given by equation ([5*2]) . The "prime" effective frequency is 
defined by slightly different formula 

1 1 /■+« d£ 1 



O' / ^TW^T 7q[2^ q + n(ie,q)+n(-Cq)]+2n(^,q) ). (59) 

For small momentum 7 q «l. Therefore fl' ~ fi^ q for q <C 1. According to computations 
the relation I/O! ~ 7 q /f2^ q is valid with reasonable accuracy for arbitrary q. Let us note 
that due to the spin-liquid equation (|5T| ) /(0) = 1. Together with (|57| ) this gives correct 
value of spin: (S 2 ) = 3/4. 

The results of numerical calculation of the functions / and g are given in Fig. 12. Distance 
between sites attains the values r = \Jm 2 + n 2 . We plot lng(r) for odd value of m + n, and 
ln/(r) for even value of m + n. One can see that for r > 2 the dependence is practically 
linear. Hence 

(S(*,r) • S(t,0)> ex (-l) m+n exp(-r/e A /), /or r > 2. (60) 

The values of the magnetic correlation length £m found from Fig. 12 are given in the Table 
II. 

Using eqs.fl5"T|) and Q58D one can express static magnetic formfactor in terms of the effective 
frequencies and Q' uk 



SW(q) = £exp(iq-r)(S(t,r) • S(t,0) = -\ +2^^ ^ ) (61) 



For q = the integrand in (gl]) vanishes at small k because Vt^ = Q' vk at k -C 1. There- 
fore only large momenta contribute into 5a/(0). For q = Q = (tt, 7r) situation is different 
because fi^k+Q = ^k and fi^k+Q = — ^k- Therefore contribution of small k into Sm(Q) is 
stressed like one over spin-wave pseudo-gap squared (1/V 2 ). The values of 5m(0) and Sm(Q) 
for different hole concentrations are given in the Table II. The static formfactor Sjw(q) at 
hole concentration 5 = 0.15 is plotted at Fig. 13. We have to note that due to the strong 
dependence on the pseudo-gap which is nearly zero at 5 = 0.05, the value of SW(Q) at 
this concentration is sensitive to small variation of parameters. Therefore the accuracy of 
calculation of Sm(Q) a t this point is very poor. This is reflection of the closeness to the 
critical concentration 6 C at which system drops down into the irregular spin-glass state (see 
discussion in the previous Section). 

Although we do not observe the shift of the Sat (q) maximum from the point Q = (tt, tt) 
our result agrees qualitatively with that of finite cluster numerical simulation (see also 
review paper and references therein, where results for other models are presented). 
However the values of Sm(Q) obtained by us are substantially higher than those from finite 
cluster simulations. We have mentioned already that SW(Q) is very sensitive to the spin- 
wave pseudo-gap at small momenta. Therefore the possible reason of disagreement is that 
with relatively small claster it is impossible to reproduce correct spin-wave spectrum at small 
momenta. 



VII. CONCLUSION 



In the present work using the picture of spin-wave condensation we have derived explicitly 
spin liquid ground state of the doped two-dimensional t — J model. Translation invariant 
spin-liquid solution exists at doping 5 > 5 C ~ 3 — 4%. We have demonstrated that spin-wave 
condensation and condensation of Cooper pairs (superconducting pairing of mobile holes) 
are intrinsically connected. The developed approach allows to calculate any property of the 
ground state as well as spectrum of excitations. In the present work we have calculated 
spin- wave pseudo-gap A M (Table II), effective spectrum (or pseudo-spectrum) of S = 1 
excitations (Fig. 11), magnetic correlation length £ M and static magnetic formfactor Sj^q) 
(Table II and Fig. 13). 



VIII. ACKNOWLEDGMENTS 



We have derived the basic results presented in the section III together with M. Yu. 
Kuchiev. I am also very grateful to him for numerous helpful discussions. I thank V. V. 
Flambaum, G. F. Gribakin and A. V. Dotsenko for critical and stimulating discussions. 



REFERENCES 



a Also at the Budker Institute of Nuclear Physics, 630090 Novosibirsk, Russia 
[1] V. V. Flambaum, M. Yu. Kuchiev, O. P. Sushkov. Physica C 227, 267 (1994). 
[2] V. I. Belinicher, A. L. Chernyshev, A. V. Dotsenko, O. P. Sushkov. Phys. Rev. B, 51, 
6076 (1995). 

[3] J. R. Schrieffer, X. C. Wen, and S. C. Zhang, Phys. Rev. Lett 60, 944 (1988). 
[4] N. E. Bickers, D. J. Scalapino, S. R. White, Phys. Rev. Lett. 62, 961 (1989). 
[5] A. Kampf and J. R. Schrieffer, Phys. Rev. B 41, 6399 (1990). 
[6] N. Bulut and D. J. Scalapino, Phys. Rev. B 45, 2371 (1992). 

[7] P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. Lett. 67, 3448 (1991); Phys. 
Rev. B 46, 14803 (1992). 

[8] P. Monthoux and D. Pines, Phys. Rev. Lett. 69, 961 (1992); Phys. Rev. B 47, 6097 
(1993); Phys. Rev. B 49, 4277 (1994). 

[9] E. Dagotto, A. Nazarenko, and A. Moreo, Phys. Rev. Lett. 74, 310 (1995). 
[10] P. W. Anderson, Science 235, 1196 (1987). 
[11] B. Shraiman and E. Siggia, Phys. Rev. Lett. 62,1564 (1989). 
[12] T. Dombre, J. Physique 51, 847 (1990). 
[13] A. Singh and Z. Tesanovic, Phys. Rev. B 41, 614 (1990). 
[14] R. Eder, Phys. Rev. B 43, 10706 (1991). 
[15] J. Igarashi and P. Fulde, Phys. Rev. B 45, 10419 (1992). 
[16] O. P. Sushkov and V. V. Flambaum, Physica C 206, 269 (1993). 
[17] A. V. Chubukov and S. Sachdev, Phys. Rev. Lett. 71, 169 (1993). 
[18] A. Sokol and D. Pines, Phys. Rev. Lett. 71, 2813 (1993). 
[19] S. Sachdev, Phys. Rev. B 49, 6770 (1994). 

[20] A. V. Chubukov, S. Sachdev and J. Ye, Phys. Rev. B 49, 11919 (1994). 

[21] J. E. Hirch, Phys. Rev. Lett. 54, 1317 (1985). 

[22] C. Gross, R. Joynt, and T. M. Rice, Phys. Rev. B 36, 381 (1987). 

[23] J. Oitmaa and D. D. Betts, Can. J. Phys. 56, 897 (1971). 

[24] D. A. Huse, Phys. Rev. B 37, 2380 (1988). 

[25] O. P. Sushkov, Phys. Rev. B 49, 1250 (1994). 

[26] A. L. Chernyshev, A. V. Dotsenko, and O. P. Sushkov, Phys. Rev. B 49, 6197 (1994). 

[27] M. Yu. Kuchiev and O. P. Sushkov, Physica C 218, 197 (1993). 

[28] E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

[29] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994) and references therein. 

[30] O. P. Sushkov, Solid State Communications 83, 303 (1992). 

[31] G. Martinez and P. Horsch, Phys. Rev. B 44, 317 (1991). 

[32] T. Giamarchi and C. Lhuillier, Phys.Rev. B 47, 2775 (1993). 

[33] H. Eskes, G. A. Sawatzky, and L. F. Feiner, Physica C 160, 424 (1989). 

[34] V. V. Flambaum and O. P. Sushkov, Physica C 175, 347 (1991). 

[35] V. I. Belinicher and A. L. Chernyshev, Phys. Rev. B 47, 390 (1993). 

[36] Z. Liu and E. Manousakis, Phys. Rev. B 45, 2425 (1992). 

[37] C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 6880 (1989). 

[38] M. Takahashi Phys. Rev. B 40, 2494 (1989). 

[39] A. V. Dotsenko and O. P. Sushkov, Phys. Rev. B 50, 13821 (1994). 

[40] J. Igarashi and P. Fulde, Phys. Rev. B 45, 12357 (1992). 

[41] G. Khaliullin and P. Horsch, Phys. Rev. B 47, 463 (1993). 

[42] K. W. Becker and U. Muschelknautz, Phys. Rev. B 48, 13826 (1993). 

[43] D. Murray and O. P. Sushkov, to be published. 



[44] Sometimes we use term "hole" in sense oft — J model and sometimes (usually in combi- 
nation particle-hole) in sense of Green's function technique. In each particular case the 
meaning of term is evident from context and can not cause misunderstanding. 

[45] J. Tranquada et al, Phys. Rev. B 46, 5561 (1992). 

[46] J. Rossat-Mignod et al, Physica B 192, 109 (1993). 

[47] M. Sato et al, J. Phys. Soc. Jpn. 62, 263 (1993). 

[48] L. P. Regnault et al. Physica C 235-240, 59 (1994). 

[49] N. M. Plakida, P. Horsch, A. Liechtenstein, and V. S. Oudovenko, Preprint JINR, E- 

17-95-287, Dubna, 1995. 
[50] F. J. Dyson, Phys. Rev. 102, 1217,1230 (1956); S. V. Maleev, Zh. Eksp. Teor. Fiz. 30, 

1010 (1957) [Sov. Phys. JETP 6, 776 (1958). 
[51] A. Moreo, E. Dagotto, T. Jolicoeur, and J. Riera, Phys. Rev. B 42, 6283 (1990). 



TABLES 



TABLE I. d-wave pairing at t = 3, /3i//?2 = 7, / = 1.8 and different values of hole concentration 
5. Units correspond to J = 1. The values of l)Fermi energy €f, 2)Chemical potential at zero 
temperature /x(0), 3)Maximum value of the gap in Brillouin zone at zero temperature A max (0), 
4) Maximum value of the gap at Fermi surface at zero temperature ApmaxiO), 5) Chemical potential 
at critical temperature n(T c ), 6)Critical temperature T c . For guidance we give in brackets the value 
in degrees, assuming that J = 0.1 5eV. 
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0.097 


0.029 (50°) 


0.15 0.144 
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0.079 


0.130 


0.040 (70°) 


0.20 0.178 


0.155 
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0.097 
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0.048 (84°) 


TABLE II. Parameters relevant to spin-wave spectrum at t = 


3, ft/ft 


= 7, / = 1.8 and 


different values of hole concentration 5. Units correspond to J = 1. 


1) Renormalized particle- hole 


interaction in spin-flip channel V ++ . 


2) The parameter 


of spin- wave 


instability C defined by (41). 


3) Lagrange parameter 


v which is equal to the pseudo-gap Am in the spin-wave spectrum. In the 


brackets we give the value in degrees, 


assuming that J = 


= 0.15eV. 4) Magnetic correlation length £m 


in the units of lattice spacing. 5) Static magnetic formfactor 6*M(q) 


for q = 0. 


6) Static magnetic 


formfactor Sm{<\) for q 


= Q = (7T,7r; 


l 








5 V ++ 


c 


v = A M 


£m 


S(0) 


S M (Q) 
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0.098 (171°) 
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FIGURE CAPTIONS 

FIG. 1. Magnetic Brillouin zone and the contour plot of a single hole mean occupation 
number i£ at 8 — 0.1 and zero temperature. 

FIG. 2. Single spin-wave exchange mechanism of hole-hole attraction. 
Fig. 3. The spin-wave polarization operator in the Schrodinger representation: 
a) One spin- wave is annihilated and the other one is created. 
b,c)Two spin- waves are either annihilated or created. 

FIG. 4. Polarization operator P aa (cu,q) in one loop approximation, a - contribution of the 
normal hole Green's function, b - contribution of the anomalous hole Green's function. 
FIG. 5. Corrections to the polarization operator. a)Single spin-wave exchange, b,c)double 
spin wave exchange. 

FIG. 6. The diagram Fig. 5c represented in the Goldstone diagram technique. 

FIG. 7. Transform of the diagram Fig.6a into diagram with effective "dot" interaction. 

FIG. 8. Effective "dot". 

FIG. 9. The effective particle-hole interaction in the spin-flip channel. 

FIG. 10. Spin-wave polarization operator with chaining of effective particle hole interaction 

taken into account. 

FIG. 11. Effective spin- wave frequency Q m as a function of the bare frequency cj q . 

Solid line: the direction q oc (1,0). Dashed line: the direction q oc (1, 1). For comparison 

we present unrenormalized frequency (Q m =Uq): dashed-dotted line. 

a) : hole concentration 5 = 0.05. 

b) : 5 = 0.1. For this concentration we present also the frequency of collective excitation o q 
for the direction q oc (1,0). 

c) : 5 = 0.15. 

d) : 5 = 0.2. 

FIG. 12. In \g\ and In |/| as a functions of the distance between sites r = \Jm 2 + n 2 . If m + n 
is odd, the In \g\ is plotted. And if m + n is even, the In |/| is plotted. The curves correspond 
to hole concentrations 5 = 0.05,0.1,0.15,0.2. 

FIG. 13. Static magnetic formfactor Sm(<1) a t hole concentration 5 = 0.15. The arguments 
are q x /n and q y /n. The dashed contour corresponds to the half width. 
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